REAL IDEAL AND THE DUALITY OF SEMIDEFINITE 
PROGRAMMING FOR POLYNOMIAL OPTIMIZATION 



YOSHIYUKI SEKIGUCHI, TOMOYUKI TAKENAWA, AND HAYATO WAKI 



Abstract. We study the ideal generated by polynomials vanishing on a semialgebraic 
set and propose an algorithm to calculate the generators, which is based on some tech- 
niques of the cylindrical algebraic decomposition. By applying these, polynomial opti- 
mization problems with polynomial equality constraints can be modified equivalently so 
that the associated semidefinite programming relaxation problems have no duality gap. 
Elementary proofs for some criteria on reality of ideals are also given. 



1. Introduction 

Polynomial Optimization Problem (POP) is a problem for minimizing a polynomial 
objective function over a basic closed semialgebraic set defined by polynomial inequalities 
and equalities: 

minimize f{x) 

subject to g%{x) > 0, i = 1, . . . ,£; (1-1) 
h j( x ) = 0, j = 1, ... ,m, 

where and f,gi, hj are real polynomial functions of x 6 M n . POP represents various kinds 
of optimization problems and can be solved efficiently under moderate assumptions by 
SemiDefinite Programming (SDP) relaxations developed by several authors, in particular 
Lasserre f5] and Parrilo [11} I12j : see, for recent developments with equality constraints 
[HJ E] and references therein. 

Let M[x] denote the polynomial ring R[xi, . . . ,x n ] and be the set of polynomials 

with degree up to k. The method constructs sequences {7 k} of optimization problems and 
their dual problems {T>k} from POP (jl.ip ; 

(iPfe) minimize L(f) 

subject to L : M[x]fc — > E, linear; 

L(l) = l,L(M k ) C [0,oo) 
(Dk) maximize q 

subject to / — q S g£l, 
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where k is an integer greater than or equal to ko ■= max{[ G | gj l, \ e | ' ],deg(/)}, and 
Mfc is defined from the constraint system of POP (11.11) : 



M k 



r j h j 



i=0 



(J; e 



X 



\Tj E 



[x],deg(aigi) < k, deg(r j /i J ) < k 



Here go(x) = 1 and Sl[i] 2 is the set of the sum of square polynomials. The union 
M of all Mfc is called the quadratic module generated by g% , . . . , g% and h\ , . . . , h m . Let 
and q\ be the optimal values of and D^,, respectively. Lasserre [3] formulated 
them as SDP problems, and showed that the sequences {fl}k>k Q and {ql}k>k converge 
to the optimal value of the given POP under moderate assumptions. In addition, he 
showed that if the feasible region has nonempty interior (no equality constraints), the 
equality = holds for any k > k$; SDP has no duality gap. On the other hand, 
Marshall [7] focused on the quadratic module M from POP instead of the basic 

semialgebraic set, and then proved that ft = qt if J{K) C M holds, where K = {x E 
M n | gi{x) > 0,i = 1, . . . ,£, hj(x) = 0, j = 1, . . . , m}; the feasible region of POP CD]) and 
3(K) = {p £ M.[x] | = 0, Vx E if}; the vanishing ideal of K. However, it is difficult to 
check this assumption for a given POP. 

In this paper, we study the method through investigations on vanishing ideals of general 
semialgebraic sets in R n . Let V(I) = {x E R n \ p(x) = 0, Mp E 1} for an ideal I C R[x}. 
When we deal with a polynomial ring over C, the Hilbert's Nullstellensatz describes a 
relationship between varieties and ideals. On the other hand, for an ideal / in M.[x], the 
Real Nullstellensatz says that J(V(J)) = I if / is real; see Section 3 for the details. We give 
elementary proofs for some criteria on reality of ideals (Theorem 13.111 and Theorem 13. 12p . 
Then we discuss equivalent conditions for the equality J(5n V(/)) = /, where S is a semi- 
algebraic set (Theorem E2 and [313]). Both conditions J(V(/)) = / and J(5n V(J)) = I are 
verifiable and closely related to duality of SDP (Proposition 12. 1| ). In addition, we propose 
an algorithm to calculate generators of "J(K), using some techniques of the Cylindrical 
Algebraic Decomposition (CAD) after Collins [21 [TO]. Applying these, one can equiva- 
lently modify any POP so that associated semidefinite programming relaxation problems 
have no duality gap. No duality gap of SDP is an important property theoretically and 
practically. For example, it is one of fundamental conditions for convergence of interior 
point methods, or it is used to confirm optimality of a solution. 

This paper is organized as follows: In Section 2, we present a relationship between the 
vanishing ideal of K and the duality of SDP for POP (|1.1|) . In Section 3, elementary proofs 
are given for some criteria on reality of ideals and equivalent conditions are obtained for 
the equality 3(S Pi V(J)) = /. Algorithms for deciding reality of ideals and for calculating 
generators of J(S D V(/)) are given in Section 4. 
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2. Duality 



We discuss the duality of SDP relaxation problems for POP (II. ID . We rewrite POP 
as follows. 



where 5 = {i £ 1" 9j( x ) > 0, j = 1, . . . , and / = (hi, ■ ■ ■ , h m ) is the ideal generated 
by hi, ... , h m in R[x]. The following proposition is a sufficient condition for the duality. 

Proposition 2.1. Suppose that 'J(K) = I. Then ft = qT. Moreover if there is a feasible 
point in (Dfc), it has an optimal solution. 

Proof. If Mfc is closed in the Euclidean topology, the similar arguments in [13|. Corollary 
21] ensure ft = q^ and existence of an optimal solution. Closedness of Mj, is shown in the 
later paragraphs of this section (Theorem I2.5P . □ 

It should be noted that Marshall [7] has shown a similar result. Under the assumption 
J(-fT) C M, he showed closedness of Mk/{3{K) n Although our assumption is 

slightly stronger than his, it can be verified directly as given below. The following theorem 
gives one of verifiable conditions for our assumption. The proof is given in Section 3 and 
a decision algorithm for the condition (2) is given in Section 4. 

Theorem 2.2. Let S be a semialgebraic set in R n , I be an ideal in R[x] and let I = 
Ii Pi • • • Pi Ik be the prime decomposition of I. The following conditions are equivalent. 



(1) 3(SnV(I))=I; 

(2) For anyt (1 < t < k), dim top (5 D V(J t )) = dim/ t holds, where dim top (5 n V(J t )) 



is the topological dimension of S P\V(It). 

Here the dimension dim/ of an ideal I is defined in Section 3 (dimR[x] = dim top = 
— 1). Theorem 12.21 becomes simpler if 5° n V(I) is nonempty, where 5° is the interior of 
S. The following is a corollary of Theorem 3.13 in Section 3. 

Corollary 2.3. Suppose that I = (hi, . . . ,h m ) is prime. If there exists a feasible point 
Xo for POP (li.jp such that xq £ S° and the rank of the Jacobian matrix g(^'"''^ m ) (^o) ^ s 
equal to n — dim/, then f£ = q^. 

Example 2.1. We consider the following POP. 



where f,gi G M[x], aj G M n and bj G R. We assume that 5° n V(I) is nonempty. Then it 
follows from Corollary 12.31 that SDP relaxation problems for POP <\2.2h have no duality 



minimize f(x) subject to x G K := S n V(J) 



(2.1) 



minimize f(x) 

subject to gi(x) > 0, i = 1, . . . , t; 

hj(x) := ajx — bj = 0, j 



= l,...,m, 



(2.2) 



gap. Indeed, it is clear that (hi 
dimkerfai, • • • , a m ] = n — dim/. 



. . . , h m ) is prime. In addition, rank Jy 




= n — 
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To show closedness of we start with the following technical lemma. For d E N, 
let A(d) = {a £ N" |a| = a\ + ■ ■ ■ + a n < d}. For k > k , let di = max{d E N | 
2d + cleg ft < k}, e { = k - degh { , and X = Hf =0 K[x]^. (di) x njli^k- We define a 
mapping : X — > by 

Q = ((<70a)agA(d )' (qea)a£A(d e ),n, ■ ■ • , r m ) H-> ^^(<?ia) 2 ft + ^ r jhj- 

i a j 

It is known that ip is surjective (the Gram matrix description of sums of squares); see for 
instance [13]. 

Lemma 2.4. Under the same assumption of Proposition \2. 11 ^ ^ Qfid only if q^ a 
belongs to the quotient of ideals (I: (ft)) = {s E | sft E /} for all a E A(di) and 
i = 0,l,...,£. 

Proof. Suppose ip(q) E /. Then ^,^ a (fe) 2 ffi belongs to / and hence vanishes on K. 
Since ft(x) > 0, we have each q\ a g% = on K and thus qi a gi = 0. Therefore the assumption 
implies qi a gi E /. The converse is obvious. □ 

Let M[x]fe be endowed with Euclidean topology. The following theorem is a slight 
modification of Theorem 3.1 of Marshall [7]. 

Theorem 2.5. Under the same assumption of Proposition [2~J\ M k is closed. 

Proof. Let = {s E B^a;]^ | sft E /} and I ej = Jfl M[x] e . . Then and I e . are closed 
subspaces of vector spaces and K[x] ej . respectively. We define X = Y\ { Q Xi a xj^j -^i> 

where Xj Q = J di for a E A(c^), i = 0,l,...,£ and X,- = M[x] ej ./7 ej for j = 1,... ,m. 

Then J is a normed space. 

Let (p: X — ► Mfc/Jjt be the induced mapping by </?. Then we have is surjective, 
<^(Ag) = A 2 (^(<?) for A E R, and = {°} by Lemma EH Hence we have M k /I k = 

{Xv | A E [0,oo),v E V}, where V is the image of the unit sphere in X under (p. In 
addition, V is compact and does not contain zero element. Now we suppose p s E M k 
and p s — ► p in Let p s and p be the cosets of p s and p respectively. Then there 

exist A s > and v s E V such that p s = X s v s . By compactness of V, we may assume v s 
converges to some element v E V. Then the limit of \ s exists, since A s = ||A s v a ||/||v s || 
converges to ||p||/|M|. Therefore we have p = lim s X s v s = \\p\\/\\v\\v E M k jl k and hence 
P E M fc . 

□ 

3. Real ideals and semialgebraic sets 

The first main result of this section is Theorem I3.11| which is called the simple point 
criterion for reality of ideals and has already been provecQ by M. Marshall [8). In this 
section we give an elementary proof of this theorem and introduce another equivalent 
condition (Theorem I3.12|) . Using these, we also obtain some conditions equivalent to 

1 This assertion is also written as Proposition 3.3.16 in [l], but the proof has a gap. 
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J(V(J) n S°) = I or to 3{S n V(J)) = I for a semialgebraic set 5 (Theorem EES and 
Theorem I2.2|) . 

Let x = (x\, . . . , x n ) and = 1 or C. For an ideal / of k[x], denotes the set of all 
zero points of I in k n . For a subset V of 3(V) denotes the set of all polynomials witch 
vanish on V. A semialgebraic set in W n is a subset of the form 



where fcj G R[x] and *y is either >, > or =. 

For an ideal / of k[x], the dimension of /, dim/, is the transcendence degree of /. If 
/ is k[x] itself, dim A; [a;] is defined as —1. If I is prime, dim/ coincides with the depth 
of / (Krull dimension of k[x]/I). For an ideal / of k[x] with primary decomposition 
/ = /iD- • -H/fc, dim/ is maxj = i ... & dim v^t- The dimension does not depend on extensions 
of the coefficient field. See [8j [15] for more details. Finally, the rank of an ideal / = 
( /ij • • • > /m) is the maximal rank of the Jacobian g[f 1 1 '"''{ m j in V(J). The topological 
dimension dim top V(/) of V(/) is the maximal dimension as manifolds. 

For a polynomial f(x) = Ylafa xa C[x], / or /(x) denotes ^2 a f a x a - Note that 



Remark 3.1. The results and the proofs of this section are still valid if R and C are replaced 



by real closed field and its extension with y— T respectively. 



An ideal / in M.[x] is called real if the following equivalent conditions are satisfied. 



(2) For any integer m and any fi, ■ ■ ■ , f m G M[x], the equation X^=i ff ^ implies 
fi G / for all i. 

We give the equivalent condition for an ideal to be real. The following proposition is 
well known (see Lemma 2.5 of [4 J for example). 

Proposition 3.1. Let I = 1% n • • • fl /& 6e a primary decomposition of an ideal I C R[sc]. 
TTien / is rea/ if and only if each It is prime and real. 

For a while, we assume that ideals are prime, and investigate the reality. 

Lemma 3.2 (Theorem 1.11 of [H]). // a prime ideal I in R[x] is real, then I' = C[x]I := 
{ph | p G C[x], h G /} is prime in C[x]. 

Proof. By assumption, / C Suppose / = { fi, ■ ■ ■ , fk) C R[x] to be prime in R[x] and 

/' not to be prime in C[x]. There exist a, b,c,d<E M[x] such that a + b y^T, c + d V-T G" /' 
and 




(3.1) 



i=lj=l 



/(x) / /(x) for x G C n . 




(1) J(V(/)) = /; 



(a + 6v^T)(c + (i\^T) = J2(ui + ViVA)fi (ui,Vi G R[x]) 



and thus 



ac — bd 




Uifi and ad + be 
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These equations yield that c(a 2 + b 2 ),d(a 2 + b 2 ) are contained in /. Since c + d\f-i G" I', 
both of c and d are not in /. Hence a 2 + b 2 belongs to /. By reality of /, both a and b 
belong to /. This implies that a + by/^1 belongs to which is a contradiction. □ 

Lemma 3.3 (Lemma 3.4 of [3]). // a prime ideal I in R[x] is real, then rank /+ dim/ = n. 

Proof. Fix the generators fi, . . . , fk of /. From Lemma 13.21 /' is prime in C[x], and thus 
rank/' + dim/' = n holds, because C is algebraically closed. Suppose rank / + dim! < n 
on V(J) and let h^s be the sub-determinants of Jacobian of size r = n — dim/. From the 
assumption, there exists i such that hi is identically zero on V(I) but is not identically 
zero on V(I'). Such hi does not belong to thus hi G" /, which contradicts reality of 
/. □ 

The following lemma can be shown similarly to the proof of Lemma I3.2| 

Lemma 3.4. If I is a prime ideal in R[x] and I' = C[x]I is not prime in C[x]. Then, 
there exists an irreducible polynomial g(x) in C[x] \ I' such that g(x)g~(x) G /. 

Proof. By assumption, / C M.[x]. Suppose / = { fi, ■ ■ ■ , fk) C R[x] to be prime and I' 
not to be prime. Similar to the proof of Lemma 13.21 there exist a,b G R[x] such that 
a + b yfl V and a 2 + b 2 G I. Set g = a + b V^T. Then g G C[x] \ V and gg = a 2 +b 2 G I. 

If <7 is factorized as g = g±g2, then <7i,<?2 /' and gg = (gigi) (5252 ) G /. Since / is 
prime, g{gx or 5252 belongs to /. Reset g as suitable one. □ 

Example 3.1. (i) The ideal J = ((x 2 + y 2 )z) is decomposed as J = ( a; 2 + y 2 ) n ( z ). 
Here, / = ( x 2 + y 2 ) is prime, while V is decomposed as V = ( x + y \f-\ ) n ( x — y \/^T ) , 
thus / is not real, neither is J. 

(ii) The ideal I = ( x 2 + y 2 + z 2 ) is prime and I' is also prime. However, since the rank 
of I is zero, / is not real. 

(hi) The ideal / = (x 2 + y 2 ,z 2 + w 2 ,xz + yw,xw — yz) is prime in M[x], while /' is 
decomposed as /' = ( x + y V-T, z + w \f-\ ) n ( x — y z — w y/A. ) , thus I is not real. 
This example is essential for the proof of the following proposition. 

The following proposition is essentially a special case of Theorem 9.3 of [S], which is 
related to going-up theorem. We give an elementary and constitutive proof. 

Proposition 3.5. If I is a prime ideal in R[x] and I' = C[x] is not prime in C[x]. Then, 
there exist irreducible polynomials gi(x), . . . ,gk{x) G C[x] \ I' such that g\{x)~g~i(x), . . . , 
9k(x)gT(x) belong to I and V = ( I',gi, . . . ,g k ) n(I',gi, . . . ,~gT) (denoted as I' g n IL) is the 
prime decomposition of I' . 

Proof. From Lemma 13.41 there exists an irreducible polynomial g\ G C[x] \ I' s.t. 
gi(x)g~i(x) G I and I' = (I',gx)n(I',gi). If we set g\ = a\ + b\ V-l, since / is prime, we 
have ai, 61 G R[x] \ /. 

We decompose I' inductively with respect to k. 
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Suppose /' is decomposed as 

i' = i' g niL (3.2) 

(/' C V g , IL) as the assertion. If we set gi = ai + b{ V-i, since / is prime, ai,b{ G R[x] \ 7. 
If both and Ij are prime, the assertion follows. Suppose I' g is not prime (which implies 
l!g is not prime either). There exist g k +i, h k +i G C[x] \ Ig s.t. g k +ih k +i G We have 

g k+ ih k+1 g k+1 h k+ i £ I g Pi 4(1%] = J, 

5fc+i5fc+T G R[x] and h k+1 h k+1 G 

Since 7 is prime, 5fc+ifffc+i or h k+ \h k+ \ belongs to I. We can assume g k +ig k +i G I. 
Similar to the proof of Lemma 13.41 we can also assume g k+ i is irreducible. Set g k +i = 
fflfc+i + V-T (o-k+i, ^fc+i G \ I). Then we have, for 1 < i < k, 

a 2 + b 2 ,a 2 k+1 + b 2 k+1 G I, 
a i a k+i ~ b 2 b k+1 ,a 2 b k+1 — a k+1 b 2 G I, 
ajflfc+i + bib k+ i or ^^+1 - bib k+ i G J, 
Ojft/c+i + a k+ ibi or aA+i - ak+ih G J. 

When aiafe + i+6i6fc + i,ai6jfc + i+Ofe + ifei G I for some i, we have (aia k+ i + bib k+1 ) 2 - (aib k+1 + 
a fc +i^) 2 = (a? " a 2 k+1 ){bj - b 2 +1 ) G /. Thus, (a? - 6?) or (a£ +1 - is in /. Since / is 

prime, the former implies a«, 6j G I and the latter implies a k +i, b k+ i G I, each of which is 
a contradiction. Similarly, aia k+ \ — bib k+ \, aib k+ \ — a k+ \bi G I yields a contradiction. 

Suppose ajOfc+i + hb k +i, aj&fc+i — G I for all 1 < £ < fc. J' is decomposed as 

I' = J H J* where 

J = ( /', at + bi v^T, . . . , a k + 6fe V-I", a k+1 + v^T), 

n ( I 1 , ax - bi V-i, ... ,a k -b k V-i, a k+ i - b k+ i V-i), 
J* = { I', ai + bi V-i, ...,a k +b k V-i, a k+1 - b k+1 V-i), 
n ( I 1 , ai - bi V-i, ... ,a k -b k V-i, a k+ i + b k+ i V-i ) . 

We show I' = J, i.e. 

(ai + bi V^l)(a» - bi V-i), (a { + bi V^i)(a k +i - b k+1 V-i), 
(ak+i + b k+ i V-i)(ai - bi V-i), (a k +i + b k+1 V-i)(a k +i - h+i V-i) 
belong to I'. This means 

a i + bf, (aifflfc+i + bib k+1 ) - (oj6fc + i - a k+1 bi) V-i, 
ajOfc+i + + (a,ib k+ i - a k+ ibi) V-i, a k+1 + b k+1 

belong to which is obvious by the assumption. Similarly, the assumption that aia k+ \ — 
bib k +i, aib k+ i + a k+ ibi G I for all 1 < i < k yields I' = J*. Exchanging g k +i and g k +i, we 
have 

/; = (I' a ,9k+i)r)(lj,g k - + j) . 
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Since I' C ( I',gi } C ( J', g±,g2 } C . . . is an ascending chain, this procedure terminates 
in finite steps. Finally, we show that if I' is primarily decomposed as (|3.2p . this is the 
prime decomposition. Suppose g m G /^ for some m > 1 and g $ I' g . Then, g m 'g m G 
IgDl^f] M[x] = I. Since / is prime, gg~ G /. By construction, we have g G I' g . □ 

Corollary 3.6. Let I be a prime ideal in M[x], then dim/' + rank/' = n holds for I' = 
C[x]I. Moreover if I' is not prime in C[x], then dim/' = diml' g = dim// and rank/' = 
rank I' g = rank // hold. 

Proof. If /' is prime, the assertion is obvious. We assume /' not to be prime. We show 
the latter assertion, which leads the former immediately. By the general theory, dim/' = 
max{dim/g, dim//} holds. Now, by symmetry, it is clear that diml' g = dim//. Next, we 
show the equation for the rank. We can assume V(I') 7^ 0. By the symmetry of 

V(l') = V(I' g )uV(lL), 

it implies V(I' g ) / 0. Set g i = ai + h VTT. If a { G I' g , then a, G I' g n //fl R[x] = I, which is 
a contradiction. Since I' g C (I g ,a,i) and I' g is prime, we have diml' g > dim (/', a$ ) for all 
i, and hence 

dim(V(/p n (V ( at ) U ■ • ■ U V ( a k ))) < dim V(l' g ). 
Using d x . (af + b})\ gi= Q = 2aid Xj (ai + h we have 

rank/' = max rank/' 

y {z|Vj,a^0} y 

= max rank/' < rank/'. 

{z|Vi,ai^0} 

The opposite inequality rank/' < rank/^ = rank// is obvious. □ 

The following lemma is very elementary, but the authors could not find it in any liter- 
ature. 

Lemma 3.7. Let k > 2 and I = I\ Pi • • • fl /& the prime decomposition of an ideal I in 
C[x]. Then the rank of I on V(/i) fl • • • fl V(/fe) is less than n — dim/. 

Proof. Without loss of generality we can assume dim/i = dim/ =: d. Suppose 
rank a . / > n — d for some x £ n^ =1 V(/j). Then, there exist fi,---,f n -d £ ^ such 
that rank x ( /1, . . . , / n _rf ) = n — d and from the implicit function theorem, the topo- 
logical dimension of V ( f%, . . . , f n -d ) m a neighborhood U x of x is d. Since /1 is prime 
in C[x], dim top V(/i) is also d, which implies that a polynomial g G I\ is identically 
zero on V ( /1, . . . , f n —d) dU x , unless otherwise dim top V(/i) < d. From the inclusion 
V(/j) Pi U x C V ( fi, . . . , fn-d ) nL/r, g is also identically zero on V(/j) n U x for all 1 < i < k. 
Since I{ is prime, g = on V(/j), which implies 5 G /«. Thus, g G n^ =1 /j = /, which 
implies I = I\\ contradicts k>2. □ 

Proposition 3.8. If L is a prime ideal in M[x\ and /' = C[x]/ is not prime in C[x]. Then 
dim / + rank I < n. 
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Proof. By Proposition 13.51 I' is decomposed as (I3.2p . By the above lemma, it is enough to 
show V(I) C V(I' g )nV(lL). Since gifR E I for all i, x E V(J) C K n implies ffi (x) = £?(x) = 
and hence x € V(ijJ) n V(l£). □ 

Lemma 3.9 (Lemma 3.9 of jBJ). For a prime ideal I in k[x], I is real if and only if 
rank J(V(/)) = rank/. 

Lemma 3.10. Assume that I is prime in R[x] and I' is prime in C[x], then dim / + 
rank / = n implies dim top V(/) = dim/. 

Proof. Set d = dim/, then dim op V(I') is d. Let x° be a point in V(/) such that rank x o / = 
n — d. Let fi,..., f n -d £ / satisfy rank ( /i, . . . , ) = n — d. By a suitable reordering of 
the variables, the equations fi = can be solved for the first n — r variables as functions of 
the last <i variables in a neighborhood of or. Let u\, . . . ,u n - r be such solution functions. 
We write x = (x n _ d+1 , . . . , x n ), u = (u\ , . . . , it n _ d ) and f(x) = f(u(x),x) for / E R[x]. 
Then 

x = (u(x), x) 

holds for all x in a neighborhood of x°. Let / belong to /', then we have f(x) = 
for all x, unless otherwise the dimension of the manifold V(I') is less than d, which is 
a contradiction. Hence x is a coordinate of the manifold V(/') in a neighborhood of x°. 
Moreover, since fx, . . . , / n _,i E M[x], x is also a coordinate of the real manifold V(/). □ 

Theorem 3.11. Let I be a prime ideal in K[x] 

i) If I is real, then I' is prime and dim/ + rank/ = n. 

mJ //dim / + rank / = n, then I is real. 

Proof. 

i) It follows from Lemma 13.21 and Lemma 13.31 

ii) From Proposition 13.81 /' is prime. We set fi, . . . , f n -d and notations as in the proof of 
Lemma f3. 101 Suppose /(x) = on V ( /i, . . . , f n -d )■ The derivatives of /(x) are zero, so 

Q _ Of _ (9/ du m Of dx m 



thus we have 



C^X^ dXyyi dXj dXyyi 3x j 

m=l J m=n—d+l J 



O /j Tl d r\ i> r\ 

J m=l J 



and hence 

df df\ fdf df 



dxi ' ' dx n J \dxi' 1 dx 



^n—d 



10... 
1... 

0... 1 



— dui 


—dui "| 


dx x _ d+1 


8x n 


—du2 


— dU'2 


dx x _ d+1 


dx n 


-du n ^ d 


-du n _ d 




dx n 



which implies rankJCV(/)) = n — d. From Lemma |3.9| / is real. □ 
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Theorem 3.12. A prime ideal I in R[x] is real if and only if dim top V(J) = dim/. 



Proof. The "Only if part follows from Theorem 13,111 i) and Lemma 13.101 

We show the "if part. Suppose dim top V(I) = dim/. If I' = C[x] is not prime, then 
V(J) is included in V(I' g ) n V(//) as in the proof of Proposition 13.81 and hence we have 
dim top V(I) < dim/. Thus /' is prime in C[x]. We show that if a polynomial / G M[x] is 
identically zero on V(J), then / = on V(I'); which implies / G /'nM[x] = / i.e. / is real. 

Let fi, ■ ■ ■ , fk generate I and denote by Q' the field obtained by the extension of Q by 
the coefficients of f,fx,..., fk- There exists a point x in V(I) C M n whose transcendence 
degree is dim/ on Q . Such a point is a generic point of V(I') on Q', from Theorem 2 of 
Chapter IV of |15| . Thus / is identically zero on V(I'). □ 

Example 3.2. (i) The ideal / = ( xy ) is decomposed as / = ( x) D ( y ). dim top V(( x )) = 
dim(x) = 1 implies (x) is real prime, and similarly (y) is real prime. Hence / is real, 
(ii) The ideal I = {y 2 — xz,x 3 — yz) is decomposed as / = J Pi {x,y), where J = 
(y 2 -xz,x 3 -yz,x 2 y- z 2 ). For each, dim top V(J) = dim{(i 3 , t A , t 5 ) ; t G R} = dim J = 1 
implies J is real prime, and dim top V(( x,y)) = dim (x,y) = 1 implies (x,y) is real prime. 
Hence / is real. 

Now we return to the semialgebraic set S. We recall that S° is the interior of a semial- 
gebraic set S in R n . 

Theorem 3.13. Suppose that S° n V(/) is nonempty and that I = I\ fl • • • PI Ik is the 

prime decomposition of I. Then the following conditions are equivalent. 

(i) 3(S°nV(I)) = I; 

(ii) For any t (1 < t < k), there exists x l in V(/^) n S° such that rank^t It = n — dim/^; 
(Hi) For anyt (l<t<k), dim top V(I t ) Pi S° = diml t holds. 

Proof, (i) (ii) Suppose that there exists t (1 < t < k) such that rank^t It < n — dim/^ 
for any x* G V(/j) n S° . From Theorem 13.111 (i), there exists x l in V(It) such that 
rank x t + dim It = n, and hence the set {x G V(/^) ; rank x If <n — dim/j — 1} is a proper 
subvarietiy of V(/t). Hence, there exists a polynomial ft identically zero on V(It) H S° and 
not identically zero on V(I t ). Thus ft It- Set / G R.[x] as / = Ils=i fsi where f s G I s \It 
for s ^ t. Then / is identically zero on V(/) n S° and / G" /, which is a contradiction. 

(ii) (iii) Suppose that there exists x* in V(/j) n S° such that rank^t It = n — dim/j. 
From the proof of Theorem 13.111 there exists a neighborhood U x t of x* in M n such that 
dim top V(I t ) n ^ =dim/ t . 

(iii) (i) It is clear that dim top V(I t ) H S° = diml t implies dim top V(/ t ) = dim/ t . The 
assertion follows from Theorem 13.121 and Proposition 13.11 □ 

We give a proof of Theorem 12.21 

Proof of Theorem\2M (i) =>• (ii) Suppose that dim top V(/ t ) n S < dim/ t for some t. From 
Theorem 13. 121 dim top V(/j) = dimlt- Thus, V(It) Pi S 1 is included in some proper subvari- 
eties of V(/ t ). Hence, there exists a polynomial /t such that ft = on V(/ t ) n S and is 

10 



not identically zero on V(It)- Thus ft It, which yields a contradiction similarly to the 
proof of Theorem 13.131 (i) => (ii). 

(ii) => (i) Replace V(J t ) n 5° by V(J t ) n 5 in the proof of Theorem EM (hi) => (i). □ 

Remark 3.2. The condition obtained by replacing S 10 by S 1 in (ii) of Theorem 13.131 "For 
any i (1 < t < A;), there exists x in V(it) Pi S such that rank^t It = n — dim 1^." does not 
guarantee (i) of Theorem l2.2l Indeed, set S = {(x,y) ; g(x,y) = 1 — x 2 — (y — l) 2 > 0} and 
h(x, y) = y, then I = (y), the origin O is in V(/) and ranko 1 = 1 = 2 — dim/. However 
3(V(I) HS) = 3(0) = (x,y)is not included in /. 

4. Algorithms for testing or guaranteeing the duality 

In this section we propose an algorithm to calculate generators of ideal 3(K) = 3(S n 
V(-F)). Applying it, one can obtain an equivalent problem to POP (II. ip such that the 
resulting problem satisfies Condition (ii) of Theorem l2.21 The algorithm uses a part of the 
cylindrical algebraic decomposition (CAD) after G. E. Collins (see [2j [10] and references 
therein for basic literature). 

The following algorithm is for detecting whether the condition of 3(K) = I holds or not. 
Note that if this condition holds, / itself should be real, because / C 3(V(I)) C 3(K) = /. 
We omit details of the CAD procedures, which are illustrated in the examples below. 

Let gx, . . . ,ge and h\, . . . , h m be defining polynomials of the semialgebraic set S and of 
generators of the ideal I respectively, i.e. {gi, . . . , gg] = {gij ; 1 < i < s, 1 < j < r^} in 

Algorithm 4.1. Input: S and /. 

(1) Compute the primary decomposition of /, / = I\ n • • • H Ik- If C[x]/t is not prime 
in C[x] for some t, then / is not real, otherwise, go to (2). 

(2) For each I t = (pi,... ,p s ) do: 

(a) Choose coordinates x = [x^, . . . , X{ d ), where dt = dim/ t , such that 1 
C(x)Ij. Here C(x) denotes the field extended by {x^, . . . ,Xi dt } from C. 

(b) Let P n and Q n denote the set of polynomials {p±, . . . ,p s } and {pi, . . . ,p s , gi, . . . , 
respectively. Execute the projection of CAD for the polynomial sets P n and 
Q n from R n to R 1 , where ieK*. For d t < n' < n let P n ' and Q n ' denote the 
set of irreducible factors of resulting polynomials on M. n from P n and from 
Q n respectively. Also let C n ' denote the set of cells in M n ' from Q n . 

(c) For any open cell U s G C dt , take a sample point x s £ U s C M. dt . 

(d) Lift x s to the point where some polynomial p G p d t+i - ls zerQ5 [_ e _ 
p(x s ,Xi dt+1 ) = 0. Denote x dt+1 = (x s ,Xi d+1 ) (x dt+1 can be more than one 
point). 

(e) Iterate the above step to the top level. Condition (ii) of Theorem 12.21 holds if 
and only if there exists a point x s such that the point can be lifted to a point 
x n e ]g>n b e i on gi n g to S. 

n 



Remark 4.1. One can compute the primary decomposition, by using computer algebra 
systems; e.g. Macaulay2, Singular, Risa/Asir. 



Remark 4.2. The condition 1 C(x)I is equivalent to the condition that the topological 
dimension of projection of V(I') C C n to C d is d. 

Example 4.1. Set h\ = y 2 — xz, hi = x 3 — yz, /13 = x 2 y — z 2 and g\ = 1 — (x — l) 2 — 
(y — l) 2 — [z — l) 2 . Applying Algorithm 4.1, we have the following: 

(1) I = (hi, hi, hs ) is prime and C[x]I is also prime. 

(2) The dimension of / is 1 and 1 ^ C(xi)I. The set of irreducible factors of the resultants 
and sub-resultants for / with respect to z is P 2 = {x, y, x 4 — y 3 }. That for / and g\ is 

Q 2 = {x, y, x 4 - y 3 , 1 - (x - l) 2 - (1 - y) 2 , 

x 2 (x - l) 2 + x 2 {y - l) 2 + y A - 2xy 2 ,y 2 (x - l) 2 + y 2 (y - l) 2 + x e 





2x 3 y, 



-1 





x 2 y 





-1 





1 


-2 


(x-l) + (y 





1 


-2 



1) 



2 

x y 


" 2 + (y 



} 



(x- \ r + (//- I) 

(Fig. [1]). The irreducible factors of the resultants and sub-resultants for Q 2 with respect 
to y is approximately 

{x,x-2,x- 0.522613, x - 1.39169, x - 0.714577, x - 1.74196, etc}. 

Let x = 0.6* as a sample point of the interval (0.5226,0,714577) E C 1 . By lifting it 
to xy-plane by P 2 , we obtain two points (0.6* , 0), (0.6* , 0.6s ). Further, by lifting it to 
xyz-plane by /, we obtain only one point (0.64,0.63,0.6s), which satisfies gi > 0. Thus, 
3(K) = I holds. 




Figure 1. Real roots of each polynomial in Q 2 in Example 14. 2 1 The results 
of projection to xy-plane by the CAD 
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Even if 3{K) ^ I, we can generate new equality constraints for POP 111. ip by the 
following algorithm. 

Algorithm 4.2. Input: S and I such that "J(K) / I. 

(1) (a) Compute the primary decomposition of I, I = I\ D • • • D Ik- 

(b) If If is not prime in M[x], replace It by its associated prime \JT t . 

(c) If = C[x]/f is not prime in C[x], I[ is decomposed as Proposition 13.51 Add 
~Rsgi,\mgi to I and go back to (a). 

(2) If necessary, operate an invertible linear transformation on W 1 so that 1 
C(xi, . . . , Xd t )It holds for all t, where dt denotes diml t . 

(3) Let to be one of the numbers satisfying J(V(it ) flS) / If and dimI to = d max , 
where d max = max{dim/f ; 1 < t < k, J(V(I 4 )n5') ^ It}- Such to can be calculated 
by Algorithm 4.1. 

(4) For each I t , if J(V(/t) n S) = It let /* be a polynomial in I t \ It , otherwise do for 
I t = (pi,...,p s ): 

(a) Let P n and Q n denote the set of polynomials {p\, . . . ,p s } and {p\, . . . ,p s ,gi, ■ ■ ■ , 
respectively. Execute the projection of CAD for the polynomial sets P n and 
Q n from M. n to R 1 , eliminating x n , ... ,X2 in this ordering. For d < n' < n, 
let P n and Q n denote the set of irreducible factors of resulting polynomials 
on R n from P n and from Q n . Also let C n ' denote the set of cells in W 1 from 
Q n . 

(b) Let {V s } denotes the set of cells in C n such that V s is included in V(It) n S. 
By assumption J(V(/ t ) n5) / It, the dimension of Proj^ V s (the projection 
of V s to M. dt ) is less than d t . Define f s as one of the irreducible defining 
polynomials of Proj dt V s in Q dt , and let /* be the least common multiple of 
/*'s for all s. Since 1 C(x%, . . . ,Xd t )It' for all t' such that df > dt, while 
1 G C(xi, . . . , Xd t ) { f t (x\, ■ ■ ■ , Xd t ) ), the polynomial / is not included in If 
for all t' such that d t > > dt. 

(5) Replace I by ( I\ , f 1 ) n • • • n ( Ik, f k ) (gi, ■ ■ ■ ,ge are not changed) and return to 
Algorithm 4.1. 

By the following proposition, the resulting ideal is strictly larger than I. By Noether 
property, these procedures terminate in finite steps. 

Proposition 4.1. Ifd(K) / I, then { h, f 1 ) n • • • n ( h, f k ) in Algorithm 4.2 is strictly 
larger than I. 

Proof. We show / = f\ t f I to , which implies / ^ / = I\ fl ■ • ■ flifc. Suppose f £ It , then 
since It is prime, /' belongs to It for some t. If J(V(It) H S) = It, then by the beginning 
of (4), /* does not belong to I to . Thus, we have J(V(I t ) n S) ^ I t . Hence by (4.b), /* does 
not belong to If for all t 1 satisfying d t < > (it, while we have supposed /* G It . Thus we 
have <i to < dt, which contradicts 0(V(I t ) H S) ^ It by the definition of to in (3). □ 
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Example 4.2. Set hi = (x 2 + y 2 + ^ 2 )(z - 2) and #1 = 1 - x 2 - (z - I) 2 . Applying 
Algorithm 4.2, we have the following. 

(1) I = I\C\ I2, where I\ = { x 2 + y 2 + z 2 ) and I2 = ( z — 2), is the prime decomposition 
and C[x]7i and C[x]/2 are prime. 

(2) Both dimensions of l\ and I2 are 2 and 1 C(x,y)I. The set of irreducible factors 
of the resultants and sub-resultants for I\ and g\ with respect to z is Qj = {x + 1 

1, x 2 + y 2 , 4x 2 + 4y 2 + y 4 }. The irreducible factors of the resultants and sub-resultants for 
Q\ with respect to y is {x, x — 1, x + 1}, thus we have 

cj = {(-00, -1), -1, (-1, 0), 0, (0, 1), 1, (1, 00)}. 

By lifting C\ to xy-plane by Q 2 , we obtain 

C 2 = {{x < -l},{x = -1},{-1 < x < 0},{x = 0,y < 0},{x = 0,y = 0} 
{z = 0, y > 0}, {0 < x < 1}, {x = 1}, {x > 1}}. 

Further, by lifting it to xyz-plane by P 3 = I±, we obtain only one point U 2 = (0,0,0) 
satisfying g\ > 0. We take defining polynomial f 1 of (0, 0) from Q 2 , e.g. f 1 = x 2 + y 2 . 

Similarly, from I2 and (71 we obtain Ql = Q% = {x, x+ 1, x— 1}. By lifting them to xyz- 
plane by P 3 = ^2, we obtain {x < 0, z = 2}, {x = 0, z = 2}, {x > 0, z = 2}, among which 
only J/f = {x = 0, z = 2} satisfy yi > 0. We take defining polynomial f 2 of {x = 0} C M 2 
from Q 2 , i.e. / 2 = x. 

We replace / by 

I = hnl 2 = (x 2 + y 2 + z 2 ,x 2 + y 2 )D(z-2,x) = (x 2 + y 2 ,z 2 )D(x,z-2) 

(this does not satisfy 'J(K) = I yet). From x 2 + y 2 = (x + y \/A)(x — y we should 

add x, y to I\. Thus we should replace / by 

I = ( x, y, z ) n ( x, z - 2 ) = ( x, y(z - 2), z{z - 2) ) . 

Now, if we set h\ = x,fi2 = y(z — 2), hs = z(z — 2) and g\ = 1 — x 2 — (z — l) 2 > 0, then 
3(K) = I is satisfied. 
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